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基于 改进 的 萤火虫 优化 算法 的 混合 语音 盲 分 离 
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摘 要 : 莹 火 虫 优化 算法 是 模拟 自然 界 莹 火 虫 彼此 之 间 通 过 荧光 素 进 行 信息 交流 的 特性 发 展 而 来 的 ， 在 搜索 复杂 函数 
全 局 极 值 点 方面 非常 有 效 。 针 对 传统 盲 源 分 离 优 化 算法 对 分 离 性 能 影响 较 大 的 局 限 性 ， 提 出 了 一 种 基于 改进 的 莹 火电 
优化 的 混合 语音 盲 分 离 算 法 。 新 算法 将 萤火虫 的 飞行 跨度 由 固定 取 值 变 为 由 新 构造 的 函数 自 适 应 调整 ， 在 加 快 收敛 速 
度 的 同时 ， 人 避免 算法 早熟 现象 的 发 生 。 实 验 结果 表明 ， 与 基于 梯度 、 标 准 营 火 虫 和 粒子 群 优化 的 盲 分 离 算 法 相 比 ， 新 
算法 对 混合 语音 信号 的 分 离 效 果 较 好 ， 在 收敛 速度 和 分 离 精度 方面 都 有 所 提升 。 
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Blind separation of speech mixtures based on improved glowworm swarm optimization 
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Abstract: Glowworm Swarm Optimization (GSO) is a population-based metaheuristic that a glowworm carries luciferin along 
with itself to exchange information with its companions. GSO is very effective in searching for the global extremes of a 
complex function. Aiming at the limitation of traditional optimization algorithms for Blind Source Separation (BSS) , this 
paper proposes a blind Speech separation algorithm based on an improved GSO (IGSO) . The new algorithm constructs a right 
function to adjust adaptively the flight span, so it cannot only accelerate the convergence speed but also avoid the premature 
convergence. The mixed speech separation experiments show that the proposed algorithm performs better than BSS based on 


Gradient, GSO or Particle Swarm Optimization (PSO) , and improves both convergence rate and separation accuracy. 
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群 智能 优化 算法 是 人 工 智 能 一 个 新 的 研究 分 支 ， 不 要 求 

标 函 数 和 约束 的 连续 性 与 凸 性 ， 在 复杂 问题 求解 方面 具有 独特 

盲 源 分 离 (blind source separation，BSS) 是 指 在 不 知 源 信 的 优势 ， 有 具有 成 本 低 、 速 度 快 和 和 鲁 棒 性 好 的 特点 ， 已 经 在 多 个 
号 和 传输 通道 特性 的 情况 下 ， 仅 由 观测 信号 恢复 出 源 信号 各 个 ” 领域 获得 了 应 用 ， 如 计算 数学 、 航 空 航天 、 交 通 、 计 算 机 科学 、 
独立 成 分 的 技术 ， 是 当前 信号 处 理 领域 的 研究 热点 ， 其 研究 成 子 通 信 、 电力、 材料 力学 和 生态 系统 等 是 ,为 了 解决 传统 BSS 
果 已 广泛 应 用 于 语音 识别 、 图 像 处 理 、 生 物 医学 和 机 械 故 障 检 人 化 算法 的 不 足 ， 群 智能 优化 算法 在 最 近 几 年 逐渐 被 应 用 于 
测 等 诸多 领域 1"。 一 般 说 来 ，BSS 算法 有 两 个 重要 组 成 部 分 : BSS， 如 遗传 算法 四， 粒子 群 优化 算法 4 和 蚁 群 算法 由 等 。 群 
标 函 数 和 和 优化 算法 。 目 标 函 数 用 来 规定 或 描述 恢复 信号 的 智能 优化 算法 属于 启发 式 算法 ,是 模拟 自然 界 蚂 蚁 , 鱼 , 蜜蜂 ， 


统计 独立 性 ， 而 优化 算法 的 任务 则 是 在 目标 函数 空间 寻找 到 最 蛙 和 细菌 等 生物 种 群生 存 或 砚 食 的 方式 ， 算 法 粒子 彼此 之 间 
优 解 〈 分 离 矩阵 )。BSS 算法 性 能 之 收敛 速度 和 分 离 精 度 主 要 ”通过 相互 竞争 和 合作 来 完成 寻 优 任务 。 

依赖 后 者 , 因此 如 何 选择 合适 的 优化 算法 是 BSS 技术 的 关键 挑 萤火虫 优化 算法 是 一 种 新 颖 的 群 智能 优化 算法 ， 是 受 自然 
战 。 传 统 BSS 优化 算法 基于 梯度 〈Gradient)， 除 非 初始 参数 值 ”，” 界 董 火 虫 通过 菊 光 素 进行 信息 交流 这 种 群体 行为 的 启发 而 构造 
选取 恰当 ， 和 否则 使 用 这 类 算法 求 得 的 解 不 会 太 理想 。 然 而 ， 的 ， 在 多 峰 函 数 寻 优 方 面 有 较 好 表现 。 标 准 莉 火 虫 算法 有 两 个 
于 言 假 设 ， 获 取 合理 的 参数 值 是 非常 困难 的 。 特 别 地 ， 当 目标 ”版 本 :一 是 由 印度 学 者 Krishnanand 等 人 提出 ,被 称 为 gowworm 
函数 空间 不 连续 或 不 可 微 时 ， 这 类 算法 甚至 无 法 使 用 。 Swarm optimization (GSO) [10;， 二 是 由 剑桥 学 者 Yang 提出 ， 
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被 称 为 Firefly Algorithm CFA) 02 131。 
同 , 但 在 具体 实现 方面 有 一 定 差异 ， 
的 研究 基于 GSO 展开 。 
为 了 克服 传统 优化 算法 的 局 限 性 和 发 挥 群 智能 优化 方法 的 
优势 ， 本 文 研究 利用 一 种 改进 的 莹 火 忠 优化 算法 来 解决 混合 语 
音 盲 分 离 问 题 。 


1 ” 盲 源 分 离 基 本 算法 
1.1 数学 模型 


两 种 算法 的 仿生 原理 相 
其 中 以 GSO 较为 常见 ， 下 


sw 人 (中 eR" 是 一 组 N 维 的 源 信号 ， 


设 s(?)=[s(0),s,(0),... 


A-[aaavlsRew 为 MxN 维 泥 合 和 


Tal 
到 


攻 9 


组 W 维 


)=[#(7) ,BG (7) Wr (0 eR” 是 的 观测 信号 


n(7)=[n(D),n(D),…n (中 sR" 为 噪声 ， 一 般 情况 下 不 予 考 虑 ， 


那么 观测 信号 x(〗) 由 N 维 未 知 的 独立 源 信号 s(?) 和 混合 矩阵 A 
线性 瞬时 混合 而 成 的 模型 可 以 表示 为 

x(1)=As(1)+n(?) (1) 

盲 源 分 离 的 任务 就 是 通过 某 种 方法 求 出 分 离 和 矩阵 W，, 从 而 


通过 观测 信号 恢复 出 源 信号 ， 这 个 过 程 的 数学 模型 为 
y(OD= Wx(?) O) 
其 中 ，y(D)=[y(), 交 (yw(D] sR" 被 称 为 分 离 信号 ， 为 


源 信号 的 估计 矢量 。 当 w=M 时 被 称 为 正定 盲 源 分 离 ; 当 
时 被 称 为 超 定 盲 源 分 离 ， 当 N >M 时 被 称 为 欠 定 盲 源 分 离 。 
1.2 信号 的 预 处 理 
通常 情况 下 ， 为 了 简化 算法 的 复杂 性 和 提高 计算 效率 ， 在 
混合 信号 分 离 前 一 般 需要 对 观测 信号 x() 进行 预 处 理 。 
a) 均值 化 (centering) ,也 即 是 对 观测 信号 x 进行 中 心 化 处 
里 ， 从 而 保证 是 零 均值 的 矢量 ， 最 直接 的 方法 是 x 减 去 的 其 均 
值 矢量 E[x] 。 
b) 白化 处 理 (whiting) ,也 叫 球面 化 处 理 , 通过 白化 矩阵 V 
使 得 观测 信号 球面 化 ， 使 v= Vx 去 相关 ， 也 就 是 使 Y 的 协 方差 
和 矩 阵 为 单位 矩阵 ， 即 ELv”]=I， 从 而 Elyy |]=ELWvv Wj]=I。 
对 源 信 号 预 白化 的 方法 有 多 种 , 如 主 成 分 分 量 分 析 法 (principal 
component analysis，PCA )、 奇 异 值 分 解法 (singular value 
decomposition, SVD ) 和 利用 信号 的 自 相 关 函 数 等 , 经 过 白化 处 
里 使 算法 计算 量 接近 原来 的 一 半 ， 效 率 大 大 被 提高 
c) 目标 函数 选取 。 负 炉 和 峭 度 是 信号 处 理 领 域 常见 的 信号 
统计 独立 性 度量 方法 ， 但 在 实际 应 用 中 ， 峭 度 对 野 值 (outliers) 
极其 敏感 ， 容 易 导 致 判定 结果 不 稳定 ， 因 此 负 箭 应 用 更 为 广泛 
些 。 但 是 负 糯 不 易 被 准确 计算 ， 于 是 文献 [14] 提 出 了 一 种 负 
的 蔡 代 公式 


N<M 


李 著 成 ， 等 
其 中 ，y 是 分 离 信号 y 的 第 ; 个 分 量 信号 ， 刁 是 y 的 第 7 阶 累 
积 量 。 候 如》 的 概率 分 布 对 称 ， 那 么 -0， 这 时 所 为 
70)s 0) 人) 


其 中 :就 被 称 为 3; 的 峭 度 。 
峭 度 存在 映射 关系 。 


分 析 可 见 ， 在 一 定 条 件 下 ， 负 灶 和 
上 的 表达 式 为 
(3y)=E(y)-3E(y) (5) 


2 ”萤火虫 优化 算法 


2.1 算法 原理 
首先 在 问题 解 空 间 中 随机 分 布 若干 萤火虫 个 体 ， 每 个 莉 火 
虫 个 体 代表 目标 函数 空间 内 的 一 个 解 且 在 某 一 时 刻 具 有 对 应 的 
萤 光 素 值 ， 荧 光 素 值 的 大 小 用 来 衡量 个 体 所 处 的 位 置 好 坏 ， 即 
荧光 素 值 与 解 的 适应 度 相 关 。 
每 个 草 火 虫 个 体 在 某 一 时 刻 拥有 自己 对 应 的 感知 范围 ， 而 
感知 范围 的 大 小 由 决策 半径 确定 。 萤 火 虫 个 体 在 感知 范围 内 ， 
寻找 比 自己 莹 光 素 值 大 的 所 有 个 体 组 成 邻居 ， 然 后 沿 移动 概率 
大 的 方向 朝 比 自己 萤 光 素 值 大 的 邻居 移动 ,最 后 更 新 决策 半径 ， 
不 断 重复 这 一 过 程 直到 达到 一 定 的 迭代 次 数 或 精度 要 求 ， 这 时 
所 有 的 莉 火 虫 个 体 聚 集 到 较 亮 的 董 火 虫 周围 ， 即 搜寻 到 了 目标 
函数 的 多 个 局 部 极 值 点 或 全 局 极 值 点 ， 以 此 来 达到 在 解 空 间 内 
寻 优 的 目的 。 
莹 火 虫 优化 算法 总 体 可 分 为 : 理光 素 值 更 新 、 邻 居 选 择 、 
移动 概率 计算 、 位 置 更 新 、 诀 策 半 径 更 新 五 个 步 又。 
2.2 数学 描述 
2.2.1 营 光 素 值 更 新 
萤火虫 个 体 的 荧光 素 值 越 高 ， 那 么 它 的 吸引 力 就 越 强 ， 其 
他 个 体 向 其 移动 的 概率 就 越 大 。 董 光 素 值 与 适应 度 函 数 Fimess 
相关 ， 其 数学 表达 式 为 
(1+1)=(1—p) l(t) + rrimess(x (t+1)) (6) 
其 中 : 4 表示 莹 火 忠 个 体 i 的 荧光 素 值 ，? 表示 荧光 素 值 的 
衰退 因子 ,7 表示 荧光 素 的 增强 系数 ，x ,e R” 表示 萤火虫 ;在 
实数 空间 的 位 置 ， Fimess(%) 表示 个 体 i 在 所 处 位 置 对 应 的 适 
应 度 值 。 
2.2.2 邻居 选择 
萤火虫 个 体 i 在 其 感知 范围 内 寻找 痰 光 素 值 较 高 的 其 它 萤 
火 虫 组 成 邻居 集合 WwW (领域 );， 如 式 (7) 所 示 。 
ND)={7:a, 00) (0:4(0) <000) 0) 


二 


Et 


党 


其 中 
半径 。 


必 是 萤火虫 个 体 ; 和 7 之 间 的 距离 , ”代表 个 体 的 决策 


2.2.3 移动 概率 计算 


根据 确定 的 邻居 集合 N; 计算 轮 盘 概率 上 户 ， 比 较 每 个 户 ， 
选择 大 的 作为 目标 对 象 大 


00-40) 
" 2 1(0) -4() (8) 
EN 
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2.2.4 位 置 更 新 


依据 移动 概率 ?; ， 萤 火 虫 个 体 i 进行 位 置 更 新 : 
x,(t+1)=x,(t)+ spn (Ox) (9) 
i 1 -x,( 吕 


其 中 : 出 表示 欧 氏 距离 ，spn>0 表 示 萤 火 虫 的 飞行 跨度 。 
2.2.5 决策 半径 更 新 
莹 火 虫 个 体 i 根 据 邻居 集合 N; 更 新 决策 半径 7: 


ne+D)=min {n,max {0,2(0)+ 8(n -IN, (ON (10) 


其 中 : 8 是 调节 因子 ， 
示 领 域 国 值 。 
2.3 飞行 跨度 函数 化 策略 

标准 萤火虫 算法 尽管 引入 了 对 决策 半径 的 考虑 ， 诀 策 半径 
的 大 小 与 感知 范围 内 区 光 素 值 较 高 的 个 体 数量 有 关 ， 如 果 数 量 
多 那么 决策 半径 小 ， 反 之 决策 半径 大 ， 这 样 可 以 有 效 地 避免 算 
法 陷入 局 部 最 优 ， 但 仍然 无 法 解决 飞行 跨度 取 值 固定 带 来 的 算 
法 收敛 速度 与 分 离 精度 之 间 的 了 矛盾。 所 以 ， 本 文 对 飞行 跨度 进 
行 函数 化 构造 ， 使 其 能 够 实时 自 适 应 取 值 : 在 算法 运行 初始 阶 
段 飞行 跨度 取 值 较 大 来 增强 全 局 寻 优 能 力 ， 加 速 收敛 ， 在 算法 
运行 后 期 飞行 跨度 取 值 变 小 以 此 保证 局 部 寻 优 能 力 ， 增 加 分 离 
精度 ， 飞 行 跨度 函数 表达 式 为 

SD1(I+HJ=He "+é (11) 

1>4>0 、1>w>0 均 为 正常 数 ，& 表示 飞行 跨度 冰 值 。 


3 ”基于 改进 的 萤火虫 优化 算法 的 盲 分 离 


3.1 算法 流程 图 
图 1 为 采用 改进 的 萤火虫 优化 算法 的 盲 分 离 流程 


表示 莹 火 虫 i 的 常规 感知 半径 ，n 表 


其 中 : 


号 


读 取 观测 信号 x ， 并 对 其 进 
行 中 心 化 、 白 化 处 


参数 初始 化 、 在 解 空 i 
随机 分 布 萤火虫 群 


火 虫 个 体 的 
值 4 


照 轮 盘 法 决定 草 
火 虫 个 体 的 移动 
方向 ， 在 此 基础 上 
新 其 位 置 x; 和 决 

策 半径 性 


对 


图 1 算法 流程 


3.2 ”算法 实现 步骤 
根据 图 1 得 出 算法 的 实现 步 又 : 
a) 读 取 观测 信号 x(?) ， 对 其 中 心 化 、 白 化 处 理 ; 
b) 初始 化 参数 p、 六 及 w、 和 *， 随 机 产生 并 个 分 离 矩阵 
作为 萤火虫 群 ， 在 搜索 空间 内 随机 初始 化 萤火虫 个 体 的 初始 
位 置 (%) 和 初始 决策 半径 (6) ;以 此 计算 该 萤火虫 的 初始 适 
应 度 Fimess[z(o)] (将 某 一 常数 比例 化 的 负 炳 作为 适应 度 函 
数 ); 


c) 依据 上 面 的 结果 计算 萤火虫 群 的 最 优 位 
应 度 Fitmess (x ) ; 
d) 式 (6) (7) 更 新 芝 火 虫 i 的 荧光 素 值 ;， 计算 莹 火 虫 i 的 邻 
居 集 合 w ; 根据 式 (8) 计算 莉 火 虫 ! 向 w 中 的 萤火虫 7 的 移动 
概率 马 ， 
e) 根 据 式 (9) (11) 分别 更 新 莹 火 忠 i 下 一 时 刻 的 位 置 *、 
决策 半径 坟 和 飞行 跨度 spn ; 
f) 比较 当前 萤火虫 群 的 适应 度 值 ， 得 到 最 大 值 。 若 此 值 优 
于 Fimess(x) ， 则 更 新 最 优 位 置 wx 和 最 佳 适应 度 Fimess (x ) ; 
g) 若 输出 达到 要 求 的 分 离 精度 或 程序 满足 了 结束 条 件 ， 程 
序 停止 运行 ， 输 出 zw 和 Fimess(x) ， 否 则 令 1=t+1， 转 到 步 又 
d) 继续 执行 。 

3.3 计算 机 仿真 实验 

通过 MATLAB 实验 ， 本 节 将 新 算法 与 基于 梯度 、 标 准 萤 
火 虫 和 粒子 群 优化 的 盲 源 分 离 算 法 进行 分 离 性 能 比较 。 


zw 和 最 佳 适 


3.3.1 实验 环境 及 参数 取 值 
通过 麦克 风采 集 3 路 语音 模拟 信号 ， 将 其 数字 化 后 采样 
20000 个 点 作为 源 信号 〈 图 2)， 其 直方 图 如 图 3 所 示 。 源 信号 
经 过 随机 矩阵 混合 后 产生 的 混合 信号 及 其 直方 图 分 别 如 图 4、5 
所 示 。 算 法 各 参数 取 值 情况 如 表 1 所 示 。 


表 1 算法 参数 取 值 


pb nn bh ¥ a oS 


取 值 0.4 0.6 0.08 6 5 0.03 0.06 0.03 


3 
2 10 T T T 
0 Ci 
2 i 2 
0 0.5 1 1.5 2 
x104 
0.05 
A A PA 
-0.05 : : 
0 0.5 1 1.5 2 
x 104 
0.2 ' ; 
0 RE 
02 ， ， 3 
0 0.5 1 1.5 2 
采样 点 数 104 


图 2 源 信号 


图 3 源 信号 直方 图 
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3.3.2 实验 结果 及 分 析 


orrrnnerrwer 一 一 一 a) 图 6 是 新 算法 在 更 新 30 次 时 的 分 离 信号 (直方 图 见 图 
-0.1 - + 到 到 \ 二 ZX 一 i VE ae 口 i 
0 l 


es 5 7)。 将 图 6 与 图 2 进行 对 比 ， 可 以 直观 地 看 出 源 信号 已 经 被 较 
和 x104 好 地 恢复 ， 而 此 时 其 他 几 种 算法 还 没有 收敛 ; 
Om b) 图 8 是 分 离 矩 阵 两 个 变量 W(11) 和 W(2.0 的 二 维 图 , 可 
0 05 1 15 2 以 形象 地 看 出 萤火虫 群 搜索 最 优 解 的 飞行 轨迹 ; 
| c) 图 9 是 几 种 算法 的 适应 度 曲线 比较 图 ,看 出 新 算法 比 基 
pe et 于 梯度 、 标 准 萤火虫 和 粒子 群 优化 的 盲 源 分 离 算 法 对 语音 混合 
”0 05 1 15 2 源 信号 的 分 离 性 能 要 好 ， 在 性 能 上 都 有 所 改善 : 收敛 速度 方面 

下 10 speedowiwm <speedoso<speedpso<speediso ， 精 度 方 
图 4 混合 信号 PrecGaaien < Precpso < Precoso < precioso 。 


d) 图 10 是 源 信号 和 恢复 信号 的 散 点 图 ， 从 图 中 可 以 看 出 
算法 收敛 后 ， 分 离 信号 的 顺序 和 相位 与 源 信号 相 比 有 的 发 生 了 
变化 ， 这 是 由 盲 分 离 的 不 确定 性 决定 的 [51。 

e) 算法 参数 分 析 。 表 1 的 这 些 参考 值 适用 于 大 部 分 应 用 场 
合 。 参 数 2 表示 荧光 素 值 的 衰退 因子 ， 若 取 值 为 1， 表 示 萤 光 
图 5 混合 信号 直方 图 素 值 无 记忆 。 /7 表示 适应 度 提 取 比 例 ， 反 映 了 算法 每 次 更 新 区 
光 素 值 的 增加 量 。8 也 被 称 为 萤火虫 邻 域 变 化 率 , 若 取 值 过 大 ， 


5 
AT A A 将 会 导致 邻 域 范围 仅 在 上 下 限 间 变化 ， 所 以 取 值 不 宜 过 大 。 
0 05 1 15 2 spn>0 表示 萤火虫 的 飞行 跨度 ， 不 宜 固定 取 值 ， 应 该 进行 函数 
0 x104 化 处 理 。n, 表示 邻居 冰 值 ， 若 取 值 过 小 ， 降 低 了 样本 的 多 样 性 ， 
民 T T T 
Or 则 可 能 使 某 些 优秀 个 体 被 误 抛 弃 ， 若 取 值 过 大 ， 则 降低 算法 收 
-10 : ! : 敛 速度 。4 表示 萤火虫 的 初始 奖 光 素 值 ， 可 以 任意 选取 。 为 
0 0.5 1 1.5 2 
104 董 火 虫 的 决策 半径 ， 若 取 值 过 小 ， 则 不 能 充分 感知 周围 的 莹 火 
10 Tr 虫 ; 若 取 值 过 大 ， 则 邻居 过 多 降低 了 局 部 搜索 能 力 ， 需 要 根据 
二 实际 情况 选择 。 算 法 中 种 群 规模 nw、 和 vw 没有 相对 稳定 的 参 
0 0.5 1 1.5 2 考 值 ， 应 根据 目标 函数 空间 的 特点 进行 合理 取 值 。 
采样 点 数 x104 
图 6 分 离 信 号 
图 7 分 高 信号 直方 图 0 20 40 ea 80 100 120 
图 9 适应 度 曲线 比较 
2 0 0.05 0.2 
0.8 上 0 间 0 AN 0 和 
-2 -= 06 -02 
0.6 5 0 5 -5 0 5 -5 0 5 
对 2 0 0.05 0.2 
Re 0 可 0 二 0 
-2 0.05 02 
-10 0 10 -10 0 10 10 0 10 
0.2 上 2 x10 0.05 02 
0 0 0 < 
0 0.2 0.4 0.6 0.8 1 EE -0.05 -02 
WI(1,1) -10 0 10 -10 0 10 10 0 10 


图 8 萤火虫 飞行 轨迹 10” 源 信号 与 分 离 信号 的 散 点 图 
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为 了 进一步 提高 萤火虫 优化 算法 的 全 局 搜索 性 能 ， 本 文 首 
先 对 原 算法 中 的 飞行 跨度 参数 进行 了 函数 化 处 理 ， 然 后 针对 传 
统 盲 源 分 离 优 化 算法 的 不 足 ， 提 出 了 一 种 基于 改进 的 萤火虫 优 
化 的 混合 语音 盲 分 离 算法 。 与 基于 梯度 、 标 准 萤火虫 和 粒子 群 
优化 的 盲 源 分 离 算法 的 对 比 实验 结果 表明 ， 新 算法 的 分 离 效 果 
较 好 ， 在 收敛 速度 和 分 离 精度 方面 都 有 所 提升 。 
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